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Abstract 

We study the behavior of the chemical reactions A + A — > A + 5 and A + A — > 5 + S (where the 
reactive species A and the inert species S are both assumed to be immobile) embedded on Bethe 
lattices of arbitrary coordination number z and on a two-dimensional {2D) square lattice. For the 
Bethe lattice case, exact solutions for the coverage in the A species in terms of the initial condition 
are obtained. In particular, our results hold for the important case of an infinite one-dimensional 
(ID) lattice (z = 2). The method is based on an expansion in terms of conditional probabilities 
which exploits a Markovian property of these systems. Along the same lines, an approximate 
solution for the case of a 2D square lattice is developed. The effect of dilution in a random initial 
condition is discussed in detail, both for the lattice coverage and for the spatial distribution of 
g read ants. 
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I. INTRODUCTION 



A rigorous description of the dynamics of the relevant macrovariables in reaction-diffusion 
systems requires a probabilistic multilevel approach retaining the essential features of the 
underlying many-body problem jH |^ . In this coarse-grained picture, typical macrovari- 
ables such as concentrations are no longer deterministic, but rather stochastic quantities. 
In a number of typical situations, the equations governing the dynamics of the mean con- 
centrations turn out to be identical with the classical law of mass action. In the absence of 
external asymmetries or of symmetry-breaking instabilities, the latter can be regarded as a 
mean-field (MF) law, in the sense that each part of the system is assumed to interact with 
the whole bulk at all times by means of an effective field which does not account for spatial 
effects. 

The above classical approach can be regarded as a good approximation as long as the 
characteristic time associated with the mean free path is short compared to the mean reac- 
tion time within the typical interaction radius. This is only the case if the system is well 
mixed at all times, either through external stirring or through fast internal diffusion. The 
opposite situation corresponds to the diffusion-controlled limit, where each reactant typi- 
cally explores a significant portion of space before undergoing a reactive collision, and the 
way in which reactants are distributed on a microscopic scale starts to become important 
for the determination of macrovariables such as global concentrations or, in the case of a 
lattice systems, the coverage in the different species. In such cases, classical MF approaches 
fail to describe the onset of inhomogeneous fluctuations induced by the intrinsic chemical 
noise of the system. Such fluctuations are nowadays directly observable at nanometric scales 
with the help of STM and FIM microscopy techniques j^, |5( and can be enhanced by specific 
geometric constraints and/or in low dimensions (e.g. catalytic surfaces), where external stir- 
ring is difficult and diffusion inefficient; eventually, they may give rise to non-classical effects 
such as memory of the initial condition, self-ordering phenomena, etc. 0]. Elucidating the 
role of geometry in this context is of great theoretical and practical interest in view of the 
recent progress in the development of nanoscale supports. 

Fluctuation-induced effects become even stronger in systems with immobile reactants, 
the object of the present paper. The particular systems we shall investigate here are the on- 
lattice reactions A+A — > A + S and A+A — > S+S with nearest neighbor interactions, where 
A and S denote, respectively, a site occupied by the reactive species ("occupied site") and 
the inert species ("empty site"), both assumed to be immobile. Popularly, these reactions are 
termed coalescence (CR) and annihilation reaction (AR) respectively. Various workers have 
intensively investigated the CR 0, Ej and the AR |, El Q H Q Q in the diffusion- 
controlled limit. Besides a series of applications for nucleation and aggregation systems 
[Til Il6| . the diffusion-controlled CR has also been recently used as a model for exciton 
fusion in polymers and molecular crystals 17, lH ljl, while the AR model provides a basic 



description for recombination processes and exciton annihilation 0, . In the immobile 
reactant limit, the AR model has been e.g. used to study free radical recombination on 
surfaces 21 1, cyclization reactions in polymers 22] and colloid deposition problems [2~ 



among other applications. Note the formal similarity of this model and models for dimer 
random sequential absorption (DRSA) 0, 25||. In such DRSA models, the deposition of a 
dimer on two empty sites is dual to the removal of two neighboring particles from the lattice 
upon reactions in the AR model, i.e. empty sites play the role of occupied sites and vice 
versa. There exists also a (less obvious) mapping between the CR model and a particular 
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case of random monomer filling with nearest neighbor cooperativity |26|, |27|, |28[. However, 
most studies concerning the above RSA models were performed for a fixed initial condition. 
Typically, the latter corresponds to a situation where all lattice sites are vacant, which in 
the dual picture of our model is equivalent to a fully covered lattice. In contrast, we shall 
consider here the general case in which the lattice is partially filled initially and study how 
this affects the subsequent dynamics and steady state of the system. 

Previous studies have shown that in the immobile reactant limit the one-dimensional 
CR and AR models with nearest neighbor interactions are characterized by an exponential 
decay of the mean coverage 9{t) in the reactive species to a nonergodic set of invariant 
states, as opposed to the empty state predicted by the MF equation [il||3(j. In the present 
work, we extend these results to the case of a partially filled Bethe lattice with arbitrary 
coordination number. In such loopless lattices, the relevant hierarchy of probabilities can 
be truncated exactly using a shielding (Markovian) property of the conditional probabilities 
for the state of a given site. This method is used to generalize previous results by Evans 
j3j| and by Majumdar and Privman js^]. Next, we treat the case of a 2D square lattice by 
performing an expansion based on the shortest unshielded path approximation developed 
by Nord and Evans for a series of RSA models [33[ . The results for the asymptotic coverage 
are then compared with Monte Carlo (MC) simulations. The effect of the initial condition 
on the spatial ordering induced by the reaction is also discussed by studying the dynamical 
behavior of the conditional probabilities and the associated fluctuation correlations. In the 
last section, we summarize the main conclusions and outline possible extensions of our work. 



II. THE CR AND THE AR MODEL IN BETHE LATTICES: EXACT SOLUTION 
VIA SHIELDING PROPERTY 

As a starting point, we consider an ensemble of Bethe lattices with coordination number 
z (the case z = 2 corresponds to an infinite ID lattice). In each lattice, sites are initially 
occupied at random with probability p (equal to the initial lattice coverage 9(0)). We then let 
the particles interact according to the CR (AR) scheme with nearest neighbor interactions 
specified above. By construction, the resulting statistical system will be translationally 
invariant at all times. Let us absorb the reaction rate R into the time scale by introducing 
the dimensionless time variable r = Rt. Let us denote by -Pfe(r) the probability that k 
randomly chosen nearest neighbor sites in a given lattice are all simultaneous ly o ccupied 
(fc-site cluster). The evolution equations for the ensemble probabilities Pj~ read [32J 

J^L = _ i) Pk _ "( zk _ 2k + 2) P k+1 , A; = 1,2,... (1) 

GST 2 

where v = 1, 2 for the CR and the AR model, respectively. 

The first term on the right hand side represents the destruction of a fc-site cluster by 
interaction between two particles inside the cluster and is proportional to the number of 
internal bonds (= k — 1 bonds). The second term represents the destruction of a fc-site 
cluster due to the disappearance of a particle inside the cluster upon interaction with a 
neighboring particle just outside the cluster. Such an event is only possible if a k + 1- 
site cluster preexists, implying that this term is proportional to Pk+i- Its coefficient is 
proportional to the number of bonds between the k cluster sites and external neighboring 
sites (=zk — 2k + 2). The different value of v stems from the fact that the /c-site cluster 
is only destroyed if the particle vanishes upon interaction with a filled neighbor outside the 
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cluster. In the CR case, this only happens one time out of two, since in a single event the 
particle inside the cluster has the same probability of vanishing as the neighbor particle 
outside the cluster. 

We now seek a special solution of the hierarchy JI} corresponding to our initial condition. 
One easily check that in our case -Pfc(O) = p . As it turns out, the hierarchy can be exactly 
truncated after the first two equations 

-0^2, (2a) 



dr 2 
dP^ 

dr 



P 2 -u(z-l)P 3 . (2b) 



Exact truncation is possible because of the existence of a special property for the quantity 
Qif...f{T~) = Pk+i/Pk, i-e. the conditional probability that a site is occupied given that its k 

nearest neighbors along any irreversible path starting from the site are occupied. Here " • " 
denotes an occupied site, while " i " denotes a conditioning occupied site. Rewriting Eqs. 
(J2J) in terms of these probabilities, one has 

d z 

— InQ. = -v-Q % ^ (3a) 

^Q.i = -l + yQ.r"^- 1 )^ ( 3b ) 

where the notation Q = P\ has been used. Now, in a Bethe lattice a conditioning site 
specified as occupied "shields" , i.e. clusters belonging to disjoint irreversible paths starting 
from the occupied site evolve independently of each other x . As a consequence, one has 

H 



Q*i = Q*M = Q*W- » etc - ( 4 ) 

i.e. the memory of the system is limited to the nearest neighbor site. In this sense, the 
subset of occupied sites can be said to display (first-order) spatial markovianity 35| . 

Using the shielding property (jH), Eqs. (j3J) become a closed two- variable system. Since 
the system is translationally invariant, the local probability Q is identical with the global 
coverage 6. The solution of Eqs. © reads 



Q # ( r ) = 9(r)=p 



z/(z-2) 



(5a) 



Q.+(t) = — _ v f . (5b) 



1 In the special case of a ID lattice (z = 2), this means that sites to the left of the occupied site do not 
"see" those on the right hand side. 
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Thus, the global coverage attains the asymptotic value 

-z/(z-2) 

(6) 

Regardless of the value of z, 9s increases monotonically with increasing p for the CR 
model, while in the AR model it follows a non-monotonic behavior with a universal maximum 
at p — |. This may be related to the fact that the mean asymptotic number of particles 
yielded by islands (=disconnected clusters) created by the ong oing reactions does not grow 
monotonically with the island size, as opposed to the CR case |36|, |37( . Thus, larger islands 
characteristic of high values of the initial coverage p may eventually yield a lower number of 
particles than smaller ones. In particular two-particle islands are known to disappear from 
the system, while one-particle islands survive forever. 

As expected, 9s decreases strongly as a function of z, approaching a zero value when 
z — > oo. This is in agreement with our intuitive expectation that the system must approach 
the classical MF prediction with increasing connectivity. 

Eq. (|5a|) generalizes previous results by Evans for the DRSA problem equivalent to the 
p = 1 case [31[ and Majumdar and Privman for the v = 2 case 32]. The special case of a 
ID lattice is obtained when z — > 2 + . In this limit, Eqs. © become 



9 S = 0(oo) 



P 



1 + 



z-2 



up 



9{t) = pexp(up[e T — 1] ), (7a) 
Q.+{r) = pe~\ (7b) 

The conditional probability Q a j is the same for both reaction schemes, as opposed to the 
^-dependent lattice coverage. The latter approaches the nonvanishing asymptotic value 

9 s = 9(oo)=pe- vp . (8) 

in contrast to the prediction of the MF solution 9mf{t) = p/{l + upr). Note that for v = 2 
and p = 1 the asymptotic coverage predicted by this equation is compatible with Flory's 
famous e -2 prediction for the isomorphic dimer filling problem [38L l39|. 



III. APPROXIMATE EXPANSION VS. "EXACT" MONTE CARLO RESULTS 
ON A 2D SQUARE LATTICE 

We now turn to the task of finding a suitable approximation scheme for the AR and the 
CR in lattices containing loops. The complex topology of particle clusters does not allow for 
an exact solution in this case. One must therefore resort to truncated expansions in terms 
of conditional probabilities and to MC simulations. 

As in the ID case [30], an Ursell expansion of the cluster probabilities using fluctua- 



tion correlation functions does not yield good results here either, since multisite fluctuation 
correlations do not decay monotonically with increasing number of sites and are therefore 
non- negligible (see subsection IIIIB|) . On the other hand, neglecting the cluster probabilities 
beyond a certain order may provide good agreement with simulations in the low-p regime, 
but the agreement is much worse when p « 1. In order to obtain a reasonable agreement in 
the whole range of p values, we must refine the truncation procedure inspired by a (spatial) 
Markovian property of the system analogous to the one observed for the ID system: while 
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in ID a single occupied site disconnects the reactive dynamics in the left and right half lines, 
in 2D an infinite line of sites specified as occupied decomposes the lattice into two inde- 
pendently evolving sublattices. More generally, it can be shown that for two-site processes 
such as the CR and the AR model in regular lattices of arbitrary Euclidean dimensionality 
"hyperwalls" of thickness equal to one lattice site shield one side of the lattice from the 
other El . 



Recently, Nord and Evans used this generalization of the Markovian property as a starting 
point to devise an expansion scheme in terms of conditional multisite probabilities for a 
DRSA model on a square lattice [33|. The main idea is to neglect the influence of conditioning 
occupied sites beyond a certain cut-off distance d c (measured in lattice spacings). However, 
the calculation of d c should be tailored so as to reflect the shielding property of occupied 
sites. The effective distance between an occupied site and a i site should be defined as 
the shortest "unshielded" path which is not blocked by other i sites. For instance, in the 
particular case to be studied here (an infinite 2D square lattice), the distance d c between the 
• site and the righmost i site associated with the (translationally invariant) probabilities 
Q m f, Q m ff and Q.^ is respectively 1, 4 and 5 lattice spacings. 

The case studied by Nord and Evans corresponds to the AR model with an initially 
full lattice. We shall now extend their calculations to the AR and the CR models with an 
arbitrary initial lattice coverage p. The starting point to perform the expansion are again the 
evolution equations for clusters of occupied sites. The first few equations for the evolution 
of low-order clusters are 



dP m 

dr 
dP.. 

dr 
dP. 



dr 

dP.. 
• 

dr 
etc. 



-2*/ P.., (9a) 

-P..-vP...-*"P.« (9b) 

-2I/P..-2I/P.. , (9c) 

-1P..-VP...-VP...-VP.. ~v p .* ( 9d ) 



Rewriting the first four equations of this hierarchy in terms of the Q's, we obtain 

d 



d:T 
d 



InQ. = -2uQ 9 f, (10a) 



^rlnQ.,; = -l-vQ.M-2vQ 9 J + 2vQ m j, 
jUnQ.^ = -l + vQ m u+2vQ, 



v 



(10b) 



(10c) 



(lOd) 
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A. First-order approximation 



Let us first consider the first-order approximation, i.e., we neglect those 4 sites beyond a 
distance further than one lattice spacing. We then have Q m ff—> Q 9 ^Q 9 ^ ~^ an d Eqs. 

(JTUJ) lead to the closed set of equations 

^lnQ. = -2uQ 9 f, (11a) 

^InQ.^ = -l-vQ.+ (lib) 

Taking into account the initial condition Q % (0) = <5 # ^(0) = p, these equations are readily 
integrated to obtain 

g.M = Mr)= (1 + ^_ e - r)r 

W = l + (12b) 

Notice that this result for the lattice coverage P m and the conditional probability Q m j is 
identical with the exact result in a Bethe lattice with z = 4 (cf. Eq. (J3J). 
The asymptotic result 

" s - p - (0O) = irri^ (13) 

can be expanded in powers of p to obtain 

S = p-2up 2 + 3u 2 p 3 + 0(p A ). (14) 

The different terms in the right hand side are recovered by a somewhat rougher truncation 
scheme neglecting all cluster probabilities involving more than a given number of sites k max 
in equations Q. Note that the term in p 2 in the right hand side of the formula (fTljl for the 
dilute case contains an additional factor 2 with respect to the expansion of the ID result 
(JHJ) for small p. This suggest that, in a hypercubic lattice with coordination number z, 
the pref actor of this term might have the form uz/2, as is the case in the corresponding 
expansion of Eq. (|5a|) for a Bethe lattice. 

Let us now compare the asymptotic values of the coverage 9s obtained from the first- 
order truncation of the cluster hierarchy with "exact" results from MC simulations. The 
MC algorithm for the AR and the CR model is performed as follows. At the beginning of 
each statistical realization the sites of a N x N periodic square lattice (torus) are randomly 
filled with particles until a predetermined global coverage p is attained. The elementary 
time step 5t is chosen in such a a way that each lattice site is visited once on average after 
one time unit At, i.e. St = At/N 2 . At each time step, a site i and one of its four nearest 
neighbor sites i' are chosen at random. If both are occupied, the reaction step takes place 
with probability p^ = R5t = 5t (where 5t = At/N 2 ), i.e. the particle at site i is removed 
from the lattice in the CR case, while in the AR model both sites i and i' are vacated. 
Fortunately, the convergence is rather fast with increasing linear size N and number of 
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statistical realizations n rea i. The limiting values 9s given in the tables HI and ITT1 correspond 
to N = 200 and n rea i = 5000 and its accuracy is at least equal to 10~ 5 . 

Fig. shows the coverage 9 as a function of the dimensionless time for the CR model. 
For an initially full lattice, a comparative plot between the ID solution, the MC result on 
the 2D square lattice and the simple MF approach is displayed. In the square lattice case, 
the mean coverage does not significantly change for times r « 10 and above, and its limiting 
value is found to be 0.2549, about 30.7% smaller than the ID result e -1 ~ 0.3679. As 
expected, the higher connectivity of the 2D lattice (z = 4, in contrast to z = 2 for the ID 
case) leads to an increased number of reactive events per occupied site, and the system gets 
closer to the empty state. As in the ID case, the long-time decay to the final state appears 
to be well fitted by an exponential. 

In the AR case, the simulation yields the exact value 9s = 0.0932 for an initially full 
lattice, off by about 31.1% from the exact value in ID. Fig. |2]shows the stationary coverage 

as a function of p for both the CR and the AR. The dependence is monotonic for the 
CR, whereas a maximum at p — 0.5 is observed in the AR case. As in the case of a Bethe 
lattice, this generic dependence on the initial coverage is likely to be robust in hypercubic 
lattices with arbitrary coordination number (cf. Fig. El for the AR case). 

For p = 1 and the CR [y = 1), we get 9s = 1/4 = 0.25 from Eq. ()13j) . which is smaller 
than the simulation value by 19.2% (cf. Fig. H]and Table whereas for v = 2, the formula 
(|13|) yields 9s = 1/9 = 0.1111, which is larger than the exact numerical value by 16.1% (cf. 
Fig. Eland Table ITT|) . Thus, for a sufficiently large p, the first-order truncation (Bethe lattice 
solution) underestimates the asymptotic coverage in the CR case and overestimates it in the 
AR case. 

On the other hand, for sufficiently low values of p the approximation gets better in both 
cases. Thus, for p = 0.5 the simulation value is larger than the approximated one by just 
0.8% for the CR (cf. Table H}. The fact that, for a given order of truncation, the accuracy 
increases monotonically with p in the parametric region p < 1 corresponding to a dilute 
system is by no means surprising: in the dilute limit the z = 4 Bethe lattice becomes a good 
approximation for the 2D square lattice, since "lattice animals" containing loops become 
rare. 

B. Second-order approximation 

For the second-order approximation, we take the whole set (fTUI as a starting point and 
make the approximations 

thereby retaining all unshielded paths with lengths smaller than or equal to two lattice 
spacings. With this approximation we get from Eqs. (|10|) 
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^-lnQ. = -2vQ 9 f, (15a) 



dr 
d 

dr 

d 

dr 

d_ 

d^ 



\nQ 9 j = -l + vQ 9 f-2uQ,j, (15b) 

In Q.J = , (15c) 

lnQ ^ = -2 + uQ^-2uQ^ + uQ^ . (15d) 

An analytical solution for these equations does not seem possible, but they can be in- 
tegrated numerically. The results for the stationary coverage are given in Tables |J and 
They are significantly better for the CR case; the deviation from the numerical result is 
maximal for p — 1 and is about —0.4%; its absolute value \A9s/0s\ diminishes monotoni- 
cally with decreasing p. In contrast, the maximal deviation for p — 1 in the AR case makes 
about 5.3% (cf. Fig. 

Better approximations can be obtained at higher orders, but the number of conditional 
probabilities to be taken into account grows dramatically. It then becomes necessary to 
automate the generation of the hierarchical equations. For instance, to third order one has 
24 different probabilities and to fourth order, 766 |33| . 

Nevertheless, the approximate conditional probabilities obtained from the second-order 
hierarchy ()15j) are already in good agreement with exact simulation results, both at the 
level of the stationary coverage and at the level of the time evolution (data not shown). 
Interestingly, the dynamics turns out to be qualitatively different depending on the value 
of v. In the CR case the inequality Q 9 > Q 9 j holds for all times, while in the AR case 
this is only true provided that the initial coverage is sufficiently low, i.e. for p < 1/2. This 
behavior is observed in Fig. [?l which also displays the time evolution of the other conditional 
probabilities (for typographical reasons, the symbols Qi, Q2, Q3 and Q3 used in the legend 
represent respectively the quantities Q 9 ,Q,f,Q,^ and Q j ). 

In contrast, above p = 1/2 the AR system displays a crossover between a short time regime 
for which Q # < Q 9 i and a long time regime with Q # > Q 9 j beyond a p-dependent crossover 
time (see Fig. 0). I.e., for short times the probability to find a site occupied given that its 
neighbor is occupied is larger than for a randomly chosen site with no previous information 
on the state of the neighbor site, whereas for long times the opposite is true. Remarkably 
enough, the qualitative behavior of both reaction schemes appears to be universal, in the 
sense that it remains the same in Bethe lattices of arbitrary coordination number (cf. Eqs. 

EJ). 

As far as higher order conditional probabilities are concerned, the inequality Q # > > 
Q / > holds at all times both in the CR case and in the dilute AR case with p < 1/2 

(cf. Fig. u\). However, for sufficiently short times we again observe a departure from this 
behavior at higher p in the AR case. In fact, all three conditional probabilities Q # ,Q 9 f,Q 9 A 

become larger than Q for a sufficiently large p (cf. Fig. |S]). In this regime, the detailed 
behavior of the above Q probabilities with respect to each other is rather complex and shall 
not be further discussed here. 

In order to interpret some of the above results, let us first characterize the occupation 
of a given site % by an occupation number rii (equal to one if the site is occupied and zero 
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otherwise). The fluctuation Srii is defined as the deviation from the average occupation in 
a given statistical realization, i.e. 5rii = rii — (rii). A special kind of two-site fluctuation 
correlation is then f m = ( 5rii Sn i+m ) , where i and i + m are two sites separated by m bonds 
along a ID path. By definition, f m is translationally invariant and depends only on the 
distance m. 

The behavior of the conditional probabilities in our hierarchy is given by the cluster 
probabilities P. In turn, the latter are related to the fluctuation correlations, which measure 
the reaction-induced ordering in the system. E.g., the sign of the difference Q # — Q m i is the 
same as the sign of the nearest neighbour two-site fluctuation correlation f\ = P -# — P ' 2 = 
(Q.f — Q. ) Q a ■ In all cases f% < as r — > oo, since two-particle clusters disappear. In Fig. 
E]MC computations for the dynamical behavior of the two-site correlations fi, fi-, jz and the 
three site correlation h = (5ni5rii + i Srii+i) in the dilute AR case (p < |) are shown. As in 
the CR case, one has f\ < for all times, i.e P #< > P 2 . However, as soon as p > 1/2, one 
has fi > for sufficiently short times (see FiglTUj). In other words, the probability to find a 
pair of neighboring sites simultaneously occupied is higher than if both sites are chosen at 
random. Most probably, the reason is that for short times the typical size of particle islands 
is still relatively large, and so is the value of P ; however, the react ion- induced growth of 
empty site clusters takes place at a higher rate than in the CR case. Thus, the probability 
that one finds an empty site beyond a certain correlation length from a given particle is 
comparatively high, thereby decreasing the value of P 2 . 

As for the behavior of /2,/3 and h, both schemes again display very similar qualitative 
features in the \ow-p regime. The numerical plots in Fig. El suggest that fi > and f% < 
for all times (for very short times, however, our precision does not allow to determine the 
sign of the correlation functions). In any case, this holds for the stationary values of these 
quantities as r — > oo. In terms of conditional probabilities, this means that Q % _^> Q . and 
Q.- < Q . , where "-" denotes a site in an unspecified state. Notice also that the absolute 
value of the three-site correlation \h\ becomes significantly larger than I/3I. Moreover, for 
yet smaller values of the initial coverage \h\ may get larger than |/2|. This suggests that 
any expansion of the cluster probabilities retaining only two-site correlation functions fails 
to describe the behavior, since long-range correlations propagate throughout the system in 
the course of reaction. As a matter of fact, in the ID case such an expansion leads to a zero 
stationary coverage to any order of the distance between sites |3?| . 

At higher values of p, the behavior is again modified in the AR case. The functions f'2 
and /3 change sign, and the above inequalities for the conditional probabilities change their 
direction. In contrast, h keeps its positive sign but decreases strongly. 

The analysis presented in this subsection suggests that the nature of the spatial self- 
ordering as a function of the initial condition is rather complex (specially in the AR case) 
and remains to be fully characterized. 

IV. CONCLUSIONS AND OUTLOOK 

Using the analogy with RSA problems, we have used the method of conditional probabil- 
ities to compute estimates for the lattice coverage in the framework of a unifying description 
for two different types of irreversible binary reactions, i.e. coalescence and annihilation. 
More traditional methods based on a spatial cut-off of fluctuations fail here, since the latter 
are propagated by the reactions over the whole system size. In contrast, the method of 
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conditional probabilities is exact in ID and branching media such as Bethe lattices, which 
can be used as a starting point for density expansions in other regular lattices [i]] (in the 
dilute limit, the Bethe lattice approximation should be good, since clusters with loops are 
rare). A further advantage of the method is that it provides a reasonable approximation for 
the "exact" simulation results beyond the dilute limit, thereby allowing to obtain a fairly 
good estimate in the vicinity of p = 1 (corresponding to the usual initial condition in RSA 
problems). Remarkably, the expansion for the CR model in this regime provides a better 
approximation than for the AR model. 

The approach used in the present paper can also be applied to mixed systems combining 
both coalescence and annihilation steps as well as to more complex kinds of initial condition 
[ipj . The correspondence between such models and RSA problems may prove useful in 
the context of pre-patterning of the substrate as a tool to improve self-assembly in certain 
systems. 

We have also seen that our model yields good results for the fluctuation-induced dy- 
namical behavior of the system. The main conclusion is that the subsequent dynamics of 
the spatial distribution is very sensitive with respect to the details of the initial condition, 
specially in the AR case, where several types of crossovers for the correlation functions have 
been identified. 

Possible extensions of our work include a more complete characterization of the transient 
behavior of the spatial distribution for the reactant species (and not only for the special kind 
of correlation functions considered here) as well as its dependence on the initial condition. 



However, exact decimation at any scale is in principle only possible in one dimension [41 
and probably also on Bethe lattices, but the analogous problem on a lattice with loops still 
requires the use of approximate techniques. 

In the above context, it is also of interest to compare the properties of such systems with 
those of their diffusion-controlled counterparts. This work could then be further extended 
to other systems such as the two-species annihilation A + B — > S + S. This reaction is 
known to induce reactant segre gati on at low dimensions and has been widely investigated 
in the diffusion-controlled case Hi 13 H 13 EE 

I47I l48l |. but its version with immobile 



reactants [32| has not received much attention yet. In particular, it would be interesting to 
see whether a shielding property can also be derived in this case, at least for a specific kind 
of initial conditions. 
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FIG. 1: Comparative plot displaying the analytical ID solution for the coverage, the 2D simulation 
result on a square lattice and the MF solution for the CR model. 



p 


MC Simulation 


truncation (1st. ord.) 


truncation (2nd. ord.) 


0.05 


0.04535379 


0.04535147 


0.04535233 


0.1 


0.08265606 


0.08264463 


0.08265592 


0.2 


0.1390146 


0.1388889 


0.1390138 


0.3 


0.1779778 


0.1775148 


0.1779583 


0.4 


0.2050887 


0.2040816 


0.2050757 


0.5 


0.2239603 


0.2222222 


0.2223642 


0.6 


0.2369787 


0.2343750 


0.2369514 


0.7 


0.2456631 


0.2422145 


0.2456311 


0.8 


0.2510806 


0.2469136 


0.2510572 


0.9 


0.2540033 


0.2493075 


0.2539584 


1 


0.2549411 


0.2500000 


0.2548402 



TABLE I: A comparison between the values of 9s obtained from MC simulations and the approx- 
imated values obtained by truncation of the hierarchy for the CR model. 
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FIG. 2: Comparative plot showing the p-dependence of the asymtotic coverage for the CR and 
the AR models on a 2D square lattice. 



p 


MC Simulation 


truncation (1st. ord.) 


truncation (2nd. ord.) 


0.05 


0.04131765 


0.04132231 


0.04132796 


0.1 


0.06952005 


0.06944444 


0.06950689 


0.2 


0.1025630 


0.1020408 


0.1025378 


0.3 


0.1185036 


0.1171875 


0.1184757 


0.4 


0.1255734 


0.1234568 


0.1255286 


0.5 


0.1274785 


0.1250000 


0.1274201 


0.6 


0.1259282 


0.1239669 


0.1259275 


0.7 


0.1217432 


0.1215278 


0.1219250 


0.8 


0.1150304 


0.1183432 


0.1158427 


0.9 


0.1057193 


0.1147959 


0.1078826 


1 


0.09318323 


0.1111111 


0.09812664 



TABLE II: A comparison between the values of 9s obtained from MC simulations and the ap- 
proximated values obtained by truncation of the hierarchy for the AR model. 
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FIG. 3: Nonmonotonic behavior of the final coverage as a function of the initial coverage p in 
the AR case. The continuous line displays the ID analytic result, while the dots correspond to 
simulation results on the 2D square lattice. The dot-dashed curve represents a spline fit of the MC 
results. 
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P 

FIG. 4: Exact simulation results for the CR vs. the first two orders of the shortest unshielded 
path approximation (the first order corresponds to a Bethe lattice). Note the good agreement of 
the second order results with the simulation over the whole p range. 
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FIG. 5: Exact simulation results for the AR vs. the first two orders of the shortest unshielded 
path approximation (the first order corresponds to a Bethe lattice). Note that the agreement of 
the approximate results with the simulation in the saturation region p « 1 is worse than in the CR 
case. 
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FIG. 6: Relative error of 9s as a function of p in the second order approximation for the AR case. 
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FIG. 7: Time evolution of the conditional probabilities obtained from Eqs. ()15|) for the AR case 
in the dilute regime (p = 0.25). 




FIG. 8: Time evolution of the conditional probabilities obtained from Eqs. (|15|) for the AR case 
in the saturation regime (p = 1). The legend uses the same notation as in Fig. [7| 
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FIG. 9: Dynamical behavior of the first three two-site fluctuation correlation functions and the 
three-site correlation function h in the dilute AR case (p = 0.25). For this computation we have 
performed simulations over 5000 realizations on a 200 x 200 lattice. 
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FIG. 10: Dynamical behavior of the first three two-site fluctuation correlation functions and the 
three-site correlation function h in the AR case with an initially full lattice (p = 1). The simulation 
parameters are the same as in Fig. |§1 
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